COMPUTING QUOTIENTS OF THE BRUHAT-TITS TREE FOR 
GL 2 (Q p ) BY QUATERNIONIC GROUPS 



CAMERON FRANC AND MARC MASDEU 



Abstract. In this paper we describe an algorithm for computing certain 
quatcrnionic quotients of the Bruhat-Tits tree for GL2(Q P ). These quotients 
are of arithmetic interest as they describe bad fibers of integral models of 
Shimura curves. 
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1. Introduction 

The goal of this paper is to describe an algorithm for computing quotients of 
the Bruhat-Tits tree for GL^Qp) by discrete subgroups of GL^Qp) arising from 
definite quaternion algebras over Q which arc unramificd at p. Such quotients 
are of arithmetic interest as, by work of Cerednik and Drinfeld (see [BC91 for a 
detailed account), they describe special fibers of integral models of Shimura curves 
at primes of bad reduction. 

In |Gre06[ Appendix] , Matthew Greenberg describes an algorithm for computing 
quotients of the Bruhat-Tits tree arising from subgroups of a maximal order in 
the rational Hamilton quaternions. This allows him to describe Shimura curves 
attached to indefinite quaternion algebras over Q of conductor 2p, for p an odd 
prime. Greenberg's method relies crucially on the fact that the quotients under 
consideration have exactly two vertices. As this is not always the case for other 
quaternion algebras, it is desirable to have a more general method. 

Our algorithm treats arithmetic groups arising from arbitrary Eichler orders 
within arbitrary quaternion algebras over Q. More precisely, the results in this 
paper prove the following result. 
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Theorem 1.1. Let p be a finite prime and let T denote the Bruhat-Tits tree for 
GL2(Qp). LetB/Q denote a definite quaternion algebra of discriminant N~ which 
is split at the finite prime p. Let N + be a positive integer coprime to pN~ . Let 
R C B be an Eichler X-order of level N + , and let T denote the subgroup of elements 
of reduced norm 1 in R[l/p\. Let T act on T via a splitting 

l: B®q Q p ^ M 2 (Q p ). 

Then there exists an explicit algorithm for computing the quotient T\T in time 

Q / (logff)W + 2.9) \ 

V (logp)V J ' 

where g is the genus of the Shimura curve attached to R or, equivalently, g is the 
first Betti number of the graph T\T. 

In particular, for fixed p our algorithm runs in time 

0( 5 2 (log.g) 3 ). 

We expect that the algorithm of Theorem 11.11 described in Subsection 13.21 will 
generalize to quaternion algebras over totally real fields. 

Our interest in developing this algorithm stemmed from a desire to use the quo- 
tient, and the work of Darmon-Pollack in [PS] , to compute values of rigid analytic 
modular forms and anti-cyclotomic p-adic i-functions. We will explain these and 
other applications in a forthcoming paper. 

This note is arranged as follows: in Section[2]we introduce notation and describe 
precisely which quotients of the Bruhat-Tits tree are treated in this paper. In 
Section [3] we describe the algorithm and give an estimate for its complexity. In the 
final Section|4]we show some examples of quotients. These examples were computed 
using a Sage (see |S + ll] ) implementation of the algorithm, and the code is available 
on the second author's website. 

Note that Bockle and Butenuth |BB10j have developed an algorithm for com- 
puting quaternionic quotients of the Bruhat-Tits tree for GL2(F g [T]). The authors 
wish to thank them for several helpful discussions about their work. The authors 
also wish to thank Henri Darmon, Matthew Greenberg, Victor Rotgcr and John 
Voight for their advice and encouragement. 

2. The Bruhat-Tits tree 

In this paper p denotes a fixed rational prime. 

2.1. Definition. The Bruhat-Tits tree T for GL2(Q P ) is a graph described as 
follows: the vertices of T are the homothety classes of Z p -lattices in Q 2 , where we 
regard Q 2 as a space of column vectors. We write V(T) for the set of vertices of T. 
Two vertices are joined by an unordered edge if there exist representative lattices 
Ai and A 2 for the respective vertices such that 

pA 1 C A 2 C Ai. 

If A is a lattice in Q 2 then we write [A] for the corresponding homothety class. We 
let E(T) denote the set of ordered pairs of adjacent vertices of T and we refer to 
elements of E(T) as ordered edges of T. 

Proposition 2.1. The Bruhat-Tits tree for GL2(Q P ) is a connected tree such that 
each vertex has degree p + 1 . 
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Proof. See |DT08| Section 1.3.1, Proposition 8]. □ 

The group GL 2 (Q p ) acts on 7": the action on vertices is given by multiplication 
of column vectors by matrices in GL 2 (Q p ). This preserves the adjacency relation 
described above, and so it thus describes an action on T by graph automorphisms. 
If A G GL 2 (Q p ) then we let [A] denote homothety class of the lattice in which 
is spanned by the columns of A. The map A i— [A] induces a bijection 

GL 2 (Q P )/Q p x GL 2 (Z P ) -4 V{T) 

which is GL 2 (Q p )-cquivariant for the natural left actions given by matrix multipli- 
cation. 

We would like a similar description for the oriented edges of T. For this we 
consider the oriented edge (yo , «i ) of T from the vertex vq corresponding to the 
homothety class of Z^, to the vertex 

' 1 
p 

The stabilizer of v\ for the action of GL 2 (Q p ) is precisely 



Vi = 



>'()■ 



q; 



1 
p 



GL 2 (Z P ) 



1 

p 



Thus, if wc write 



r (p) = GL 2 (Z 




GL 2 (Z P ) 



1 

p 



GGL 2 (Z P ) 



p\c 



then the stabilizer of (vq, v±) for the action of GL 2 (Q p ) is Q* • To(p). The map 

1 



A^[ [A], 



A 



p 



yields a bijection 

(2.1) GL 2 (Q P )/Q p x .r (p)4£(T) 

which is equivariant for the left action of GL 2 (Q p ). We will refer to the vertx 
as the privileged vertex and to the edge (vo,v±) as the privileged edge. 

Lemma 2.2. There is a set of coset representatives {ei}i for GL 2 (Q p )/Qp • Tq(p) 
given by matrices with coefficients in Z. Moreover, there is an effective algorithm 
that, given any matrix in g G GL 2 (Q p ), finds a corresponding scalar A G Q p and 
matrix t G To(p) satisfying: 

gXt = e 4 . 

An analogous statement holds for GL 2 (Q p )/Qp • GL 2 (Z p ). 
Proof. The matrices e\ have the form: 



V 
r 







p a 

pb r 



< r < p b+ \ 



< r < p" 
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It is easy to see what elementary matrix transformations need to be performed in 
order to reduce any matrix in A/2(Qp) to one in the above form. □ 

One can encode the vertices and edges of the Bruhat-Tits tree conveniently via 
a set of representatives as in Lemma 12.21 

Remark 2.3. It is crucial to have an efficient reduction algorithm, as in Lemma l2~2| 
if one wants to encode the T-action on T . The action of L on T is induced from 
the action on M2(Q P ) by matrix multiplication, and in order to descend this to a 
computable action on the tree, one needs to calculate the reduction of 72 for any 
7 e r and x £ M2(Q P ). This will be used in Algorithm [TJ Note also that the 
running time of such a reduction can be regarded as bounded by a constant, and 
therefore it can be ignored in the complexity estimates. 

2.2. Quaternion algebras and subgroups of interest. Let £>/Q denote a def- 
inite quaternion algebra which is unramified at p. If Hs a place of Q, then we write 
Bi = B ®q Q; where Q/ is the completion of Q at I. Thus B^ is isomorphic with 
the Hamilton quaternions and B p = M2(Qp). Let N~ denote the discriminant of 
B, so that N~ is a product of an odd number of distinct primes, and N~ is coprimc 
to p. 

Let i? max C B denote a maximal Z[l/p]-order. Let N + be a positive integer 
which is coprime to pN~ , and let R C i^ max denote an Eichler Z[l/p]-order of level 
N + . Since B is unramified at p, it satisfies the Eichler condition for Z[l/p]-orders 
and there hence exists a unique, up to conjugation by B x , Eichler Z[l/p]-ordcr of 
each level. Let t denote an isomorphism 

1: B p S M 2 (Q p ) 

which satisfies i(i?™ ax ) = M 2 (Z p ). Let L = r(p,N~,N + ) denote the subgroup of 
elements of reduced norm 1 in R. The group T acts on the Bruhat-Tits tree via 
the splitting 1. 

A discrete subgroup of GL 2 (Q P ) is said to be Schottky if it acts without fixed 
points on the vertices V(T) of the Bruhat-Tits tree. 

Proposition 2.4. (1) The group T is a finitely generated discrete subgroup of 
GL 2 (Q p ) and the quotient T\T is finite. 
(2) There exists an integer constant M > 1 depending only on pN~ such that 
if N + > M, the group T is Schottky. In this case the abelianization ofT is 
a finite free X-module of rank 

g = l-V + E, 

where V and E denote respectively the number of vertices and edges ofT\T. 

Proof. These are standard facts from the theory of Schottky groups. See, for ex- 
ample, |GvdP8Q[ Section 1.3]. □ 

Remark 2.5. Our interest in such subgroups arises from the fact that quotients 
of the form T\T describe bad special fibers of Shimura curves, cf. |BC91) . In 
the Schottky case the corresponding special fiber has irreducible components in 
bijection with the vertices of r\7". Each component is isomorphic with P 1 over F p , 
and two components meet in an ordinary double points defined over F p if and only 
if the corresponding vertices of T\T are joined by an edge. Thus, the algorithm in 
this paper describes a way to compute the bad special fibers of integral models of 
Shimura curves over Q. 
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3. Computing quotients 

3.1. Solving the T-equivalence problem. Let £>/Q denote a definite quaternion 
algebra which is unramificd at p, as in Section 12.21 Let N~ be the discriminant of 
B, so that N~ is the product of an odd number of distinct prime numbers. Let 
N + be an integer which is relatively prime to pN~ . Let R, I and T be as defined 
in Section l2~2l 

Although L is a finitely generated group, it is not defined in terms of a convenient 
finite presentation. There is, however, a natural increasing filtration of L by finite 
sets which can be used to compute a finite generating set. For each n > set 



x G R and nrd(a;) = p 



2n 



Then F„ C F n+ i since x/p n = (px)/p n+1 and 

r =U r - 



n>0 



Note that each F„ is indeed a finite set since B is definite. 

Following [BBlOj . we introduce some generic notation for group actions. If G 
is a group and X a left G-sct, and if u and v are elements of X, then we write 
Home (it, v) to denote the collection of all elements of G which map u into v. We 
also write Stabc(it) = Hoiiig(ii, u). 

The central problem that arises when computing the quotient T\T is the com- 
putation of the sets Homr (it, v) for two vertices u and v of 7" (or for two ordered 
edges) . Our solution to this problem uses short vector searches in certain Z-lattices 
of rank 4. In the following exposition of our method we concentrate on the case of 
vertices. The case of ordered edges is treated analogously. 

Assume that two vertices u and v of T are represented by reduced matrices as 
in Lemma [2.21 We will abuse notation and denote the representing matrices also 
as u and v. Define m via the formula 2m — val p (det u det v). The integer m is the 
half-length of the path that joins the vertex u to the vertex v passing through the 
privileged vertex w . Write p a = detu and p b = detw, so that 2m = a + b. Note 
that 

(3.1) Hom r (u,u) = Hom GL2(Qp) (u,u) n T 

= (wStab GL2( Q p) (u )u" 1 ) nr 

= (i;Q p x gl 2 (z p ) u - 1 ) nr. 

Lemma 3.1. If m is not an integer then Homr(u,v) = 0. Otherwise one has 

Hom r (u, v) = p- m vM 2 (Z p )u* n T, 
where u* is the matrix satisfying uu* = detu. 

Proof. Since T C SL2(Q P ), the corollary to Proposition 1 of |Ser03| Chapter 2, 
Subsection 1.2] shows that m must be an integer for u and v to be equivalent under 
r, as the group T preserves the parity of the distance between any two vertices. 
By Equation (|3.1| . it suffices to show that 

(vQ$ gl 2 (z p ) w - 1 ) n r = P - m vM 2 (z p )u* n r. 
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Write z £ (y% GL 2 (Z p )w- 1 ) n T as z = vXgu" 1 for A e Q p x and g £ GL 2 (Z P ). 
Observe that z = p~ a v\gu* . Since z £ T we have det z = 1 and hence 

\ 2 p b - a a = 1, 

where a = det g £ Z* . Therefore 2 val p (A) = a — b and thus vai p (p _a A) = — m. We 
conclude that p~ a \g belongs to p~ m GL 2 (Z p ) C p~ m A/ 2 (Zp), and thus 

(«Q£ gl 2 (z p ) m - 1 ) n r c P - m vM 2 (z p )u* n r. 

Conversely if one writes z £ p~ m vM2(Z p )u* n T as 

z = p~ m vgu* = p a ~ m vgu~ 1 , 

for g £ Ma(Zp), then one checks that det g = 1 and hence g £ GL 2 (Z p ). The 
lemma follows. □ 

Remark 3.2. Since p fixes every homothety class of lattices, it acts trivially on 
the vertices of T . Thus, the equation in the preceding lemma can be rescaled to 
read: 

Hom r (u,v) = |J p n ' m vM 2 {Z p )u* n p n T n . 

n>Q 

This expression realizes Homr(u, v) inside the Eichler order R. This is convenient, 
as although T consists of quaternions of unit norm, in practice one usually works 
with integral elements in R of norm an even power of p. 

Since Y is defined via the local splitting t, explicit knowledge of i is required to 
compute effectively with T. Of course, i is a linear map between four-dimensional 
vector spaces, and so it can be described easily as a matrix. A complicating factor 
arises from the fact that this matrix can only be stored up to a finite precision. In 
particular, the image of i is not known explicitly, which complicates the problem 
of determining if a given matrix in Af 2 (Q p ) belongs to T. In order to bypass this 
difficulty we introduce the notion of an approximation to a Q p -linear map. 

Definition 3.3. Let n > be an integer, let V and W be two finite dimensional 
Qp-vector spaces, and let Ay C V and Aw Q W be Z p -lattices. Let / : V — > W be 
a Qp-linear map satisfying /(Ay) C Aw- Then an approximation of / to precision 
?? is a Qp-lincar map g: V — > W such that g = f (mod p n ) when restricted to Ay. 

The following lemma lies at the heart of our algorithm. 

Lemma 3.4. Letu andv be matrices in Af 2 (Z)nGL 2 (Q p ) representing two vertices 
of T ■ Let f: A/ 2 (Q p ) — > B p be an approximation of i _1 to p-adic precision 2m 
relative to the orders M 2 (Z p ) and R p . Define a Z-lattice in B as follows: 

A(u, v) = f(vM 2 (Z p )u*) DR + p 2m+1 R. 

Then Homr(u,?j) is nonempty if and only if the shortest vectors in A(u,v) have 
reduced norm p 2m . 

Proof. First note that elements in A(it, v) have reduced norm of valuation at least 
2m. This is so because l transforms the reduced norm of B p into the determinant 
of A/ 2 (Qp). 

Let g £ Homr(u,v) and write 

g = p n - m vxu* = o(y), 
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where x £ M2(Z P ) and y £ p n Y n . We claim that the element A = p m ~ n f(g) is a 
shortest vector in A(u,v) of reduced norm p 2m . The first paragraph implies that 
the reduced norm has to be at least p , and since A has reduced norm p it is 
necessarily shortest. It remains to show that A belongs to A(u,v). 

Consider A' = L~ 1 {p m ~ n g) and note that A' belongs to R: clearly A' = p m ~ n y 
belongs to R[l/p\. But also A' = l~ 1 (vxu*) is p- integral, since vxu* is so, and i 
preserves p-integrality. Therefore 

A' g R[l/p] H R max p = R. 

Thus, by definition of /, we see that A' and A are congruent modulo p 2m+l and we 
can write 

A = A' + p 2m+1 a, a£R. 

This shows that A' belongs to A(u, v), and therefore so does A. 

Conversely suppose that A G A(u, v) is of reduced norm p 2m . Write A as 

A = f(vxu*) + p 2m+1 a, x £ GL 2 (Z p ), a £ R. 

Again, by definition of /, we can rewrite the above expression as 

A = r l {vxu*) + p 2m+1 a', a 1 £ R. 

We claim that t(A) belongs to Homr(u, v). First, note that i(A/p m ) is indeed an 
element of T. It remains to show that i(A) takes the vertex corresponding to u to 
that corresponding to v. For this we write 

l(X) = vxu* +p 2m+ h(a') = v(x+p 2m+1 v-h(a , )(dctuy 1 u)u*, 

and note that the matrix 

x + p 2,n+1 t)- 1 t(o , )(det u)-^ = x + p{p b v- 1 )L{a')u 

belongs to GL2(Z p ) because x does. Therefore it stabilizes the vertex vq, and this 
concludes the proof. □ 

Remark 3.5. If u and v are ordered edges of T, then the exact analogue of 
Lemma EOl is true with the Z-lattice A(u,v) replaced with 

A'(u,v) = f(vT (p)u*)r\R + p 2m+1 R. 

Remark 3.6. Lemma [3.41 is useful since the lattice A(u,v) contains the lattice 
p m+ R, and so A(u,v) can be described on a computer using p-adic approxima- 
tions for a basis, as long as the approximations have at least 2m digits of p-adic pre- 
cision. One may then use standard techniques, like the LLL algorithm as explained 
in |Coh93[ Section 2.6], to find the shortest vectors in A(u,v). These algorithms 
have a complexity of 0(m 3 ). This yields an efficient algorithm for determining 
whether vertices or edges of T are equivalent under T. 

3.2. Computing quotients of the tree. In this subsection edge always refers to 
an ordered edge. If E is a set of edges of T and if v is a vertex, we denote by E(v) 
the subset of those edges of E having origin at v. 

We have described how to determine whether edges of T are T-equivalcnt. Given 
this, it is not much more difficult to compute the quotient r\7~. The algorithm [1] 
which we describe below in fact computes a fundamental domain for the action of 
r on T, and this data is richer than the data of the quotient graph. In analogy 
with the case of Ricmann surfaces uniformized by %. a fundamental domain in T 
for the action of Y consists of: 
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(1) a connected subtree V C T whose edges form a set of distinct coset repre- 
sentatives for the action of T on the edges of T; 

(2) a collection of tuples (u, v, 7) where u and v are distinct boundary vertices 
of T>, and 7 £ T satisfies 77/ = v. We compute one such tuple for each pair 
of identified boundary vertices. The data of all the tuples is refered to as a 
boundary pairing; 

(3) in the case when T is not Schottky, a list of all the nontrivial edge and 
vertex stabilizers. 

Remark 3.7. Let T> C T denote a connected subtree whose edges form a complete 
set of distinct coset representatives for the action of T on the edges of T. If u is a 
boundary vertex of T>, then it must be T-equivalent to at least one other boundary 
vertex. If u were not equivalent to any other boundary vertices, then the edges of 
T adjacent to u and outside T> could not be T-equivalent to any edges in T>, which 
contradicts the fact that the edges of V contain a full set of coset representatives. 



Algorithm 1 Compute a fundamental domain for T acting on T 

Input: A prime p, an order R C B as above, and a splitting l p : R p = M 2 (Z p ). 
Output: A fundamental domain together with an edge pairing. 
Select a vertex v G V(T) to begin. 
Initialize a queue W with vq. 
Initialize E and P as empty lists, 
while W ^ do 
Pop v from W. 
for e 6 E(T)(v) do 

if there is no e' G E{y) which is T-cquivalent to e then 
Append e to E. 

if there is a vertex v 1 G W which is T-equivalent to t(e) then 

Append (t(e) 7 v\j) to P, where 7 G T satisfies jt(e) = v' . 
else 

Push t(e) onto W. 
end if 
end if 
end for 
end while 
return E, P. 



Remark 3.8. If one is interested in computing T-invariant harmonic cocycles on 
7", then one can facilitate this by storing extra data during the execution of the fun- 
damental domain algorithm. More precisely, to compute the harmonicity relations, 
one needs the data of how all of the edges leaving a given vertex get identified under 
r. Algorithm [T] computes this data when testing for membership in the fundamental 
domain, but with the algorithm stated as is, the data is then discarded. In prac- 
tice one can store these identifications and associate them with the corresponding 
vertices of the fundamental domain. 

Proposition 3.9. Algorithm^ terminates. 
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Proof. This follows from the finiteness and connectedness of the quotient graph 
r\7~. The algorithm incrementally explores a fundamental domain, and produces 
a boundary pairing for the boundary vertices. □ 

We wish to analyze the complexity of algorithm [T] in terms of the genus g of the 
corresponding Shimura curve. Given integer p, N~ and N + as above, denote by 
Xo(p, N~ ,N + ) the Shimura curve associated to an Eichler order of level N + in the 
indefinite quaternion algebra with discriminant pN~ . The following genus formula 
is due to Ogg, and also appears in JGR061 Prop 5.2]: 



Theorem 3.10. The genus of X (p, N~ , N+) is: 



pN~N + 
12 



n 

t\pN~ 



n 

t\N+ 



£3 

3 



£4 

4 ' 



where for k = 3,4: 



v t (k) 




l\\N+ 



n •+ (~4 



n 

l 2 \N+ 



if(T?) 

else. 



Here (7) stands for the Kronecker quadratic symbol. 
Proof. See |Ogg83[ pg. 280,301]. 

For simplicity, assume that the group T is Schottky. 



□ 



Proposition 3.11. Let n = n(p, N~ , N + ) be the number sets of the form Homr(it, v) 
that need to be computed in order to find a fundamental domain for T acting on T ■ 
Then 

(g-l)(p 3 ~2p + 2.9-1) 



n < 



(p-iy 



Proof. Let V (respectively E) denote the number of vertices (respectively edges) in 
a fundamental domain. The algorithm terminates after doing at most pE compar- 
isons among edges, and at most v ( v ~ lS > comparisons among vertices. Therefore: 



pE- 



V(V-1) 



Since we are assuming that T is Schottky, the quotient T\T is (p+ l)-regular, and 
E and V are related by the formula 



2E=(p + 1)V, 



from which we can extract: 



V = 



2(ff-l) 
p-1 ' 



E = 



(p+l)(g-l) 
p-l 



We thus obtain 



n < 



(g-l)(p 3 -2p + 2g-l) 
(p-l) 2 



□ 
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Remark 3.12. The graph T'\T (where T' C r is as in the proof above) is Ramanu- 
jan (see |LublO[ Theorem 7.3.1]) and this implies, according to |LublO| Prop 7.3.11], 
that the diameter of the graph T\T is 0{\ogg/ logp). Therefore the running time 
of algorithm Q] is 

n / (logg) 3 (g-l)(p 3 -2p+2. 9 -l) \ / (log g ) 3 g(p 3 + 2, g ) \ 

U \ (l0g P )3(p-l)2 J U { (logp)3p2 )■ 

Remark 3.13. One can reduce to the Schottky case by using the second fact in 
Proposition 12. 41 and the running time in the general case is easily seen to be 

/ (log gf{Mg - l)(p 3 - 2p + 2Mg - 1) \ = / (log Mg) 3 Mg(p 3 + 2Mg) \ 
\ (bgp)3(p-l)2 J v (logp)V J' 

where M is a positive integer depending only on pN~ , as appearing in Proposi- 
tion nu 

In practice certain improvements can be made that speed up the computation 
of a fundamental domain in the generic case where few vertices have nontrivial 
stabilizers (e.g. Schottky groups). If a vertex v of T has a trivial stabilizer in T, 
then all the edges adjacent to v are necessarily inequivalcnt under T, and so one 
need not test for membership in a fundamental domain: if one accepts the vertex 
in the domain, then one must also add all of the adjacent edges to the domain. In 
the case where T is Schottky, one can use this improvement to avoid making any 
edge comparisons in Algorithm [TJ The running time then improves to 

Q / (logff)V \ 

V(iogp)V/ ' 

Said in terms of the quotient graph, in this speed-up we are using the fact that if 
r is Schottky, then T\T is a (p + l)-regular connected multigraph without loops. 
Since for fixed pN~ all but finitely many choices of N + produce Schottky groups, 
this speed-up is significant. 

4. Examples 
We conclude with some examples of quotient graphs 

T(p,N-,N + )\T. 

As a first example we take p = 2, N~ = 13 and N + = 1. The corresponding 
graph is pictured in Figure [T] This graph corresponds to the special fiber at 2 of 
the Shimura curve of level 1 constructed from the rational quaternion algebra of 
discriminant 26. The genus of this curve is 2, as one easily deduces from the graph. 



Figure 1. r(2,13,l)\T 
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In our next example, FigureEJ we increase the level and consider T(2, 13, 9). The 
quotient map 



r(2,i3,9)\r-> r(2,i3,i)\r 



is described by the colours of the edges of T(2, 13, 9)\T. Note that this case treats 
a Schottky group, and so the graph T(2, 13, 9)\T is 3-regular. It corresponds to a 
curve of genus 13. 



Figure 2. r(2,13,9)\T 



Next we present a more complicated example, T(53, 11, 2)\T, which corresponds 
to the special fiber at 53 of a Shimura curve of genus 31. Note that due to the 
large number of edges between vertices, rather than graphing all edges we have 
simply labeled each edge in Figure |3] to indicate the number of edges that pass 
between adjacent vertices of r(53, 11, 2)\T. This example took less than 5 seconds 
to compute on an Intel Core i5 running at 3.2 GHz. 
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Finally, by way of comparison, the quotient r(211, 1511, 1)\T is a graph with 
254 vertices and 26678 edges and corresponds to a curve of genus 26425. It took 
less than 3 hours to compute on the same Intel Core i5. 
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